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Abstract 

We present results for the coupling of the light vector mesons to the tensor cur- 
rent, relative to the standard vector meson decay constants. From an (D(a)-improved 
lattice study, performed at three values of the lattice spacing in the quenched approx- 
imation, our final values (in the continuum limit), in the MS scheme at /x = 2 GeV, 
are: fj/f p = 0.72(2)(lg), /£.//*• = 0-74(2), fl/U = 0-76(1). 



1 Introduction 

1.1 Motivation for this computation 

The most direct way to extract one of the least known CKM parameters, \V u b\, is from 
studies of the corresponding leptonic and semileptonic decays of 5-mesons. Although 
the leptonic decays are still beyond the reach of the present experiments, CLEO, Belle 
and BaBar are providing rather accurate measurements of the branching ratios for the 
semileptonic B — > p{tt)Iv modes [1]. On the theoretical side, lattice QCD and light cone 
QCD sum rules (LCSR) are expected to provide model independent information about 
shapes and absolute values of the relevant form factors. That necessitates good control over 
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Figure 1: The effect of replacing fj/f p — 1 — > 0.75 in eq. (1) results in about 15 — 20% suppression of 
the form factor Af^ p (q 2 ). The ratio is evaluated by using the leading twist distribution amplitudes from 
ref. [5]. 



the low energy (non-perturbative) QCD dynamics. The complicating feature of heavy — > 
light decays is that the g 2 -region accessed by these decays is large, e.g. for B — > p£u, with 
i = e, fx, it is < q 2 < (tub — m p ) 2 = 20.3 GeV 2 . Lattice QCD results are available for 
large momentum transfers, q 2 > 10 GeV 2 (see refs. [2]), while LCSR results are expected 
to be reliable in the region of low g 2, s [3-5]. In order to reduce the uncertainties in 
the form factors obtained in the latter approach, one needs better control over the non- 
perturbative parameters which are explicitly present in the LCSR, such as moments of the 
light cone wave functions, various hadronic couplings etc. A step in this direction is made 
in this letter, where we report results for the ratio of the coupling constants of the light 
vector mesons to the tensor (fy) and to the vector current (fv), which is an important 
ingredient in the LCSR expressions for the B — > p and B — > K* semileptonic form factors. 
While the constant fy can in principle be extracted from experiments on e + e~ V° and 
r~ — > V~v T , the coupling fy can be estimated only theoretically. To do so, we will use the 
lattice QCD. The first computation of this coupling has been attempted in ref. [6]. The 
QCDSF collaboration has also presented their preliminary results in ref. [7]. Here we will 
restrict our attention to the ratio of the couplings fy/ fv, use the data generated at three 
values of the lattice spacing and extrapolate to the continuum limit. Before we enter into 
details of the lattice computations, let us make a brief assessment of the importance of 
having an accurate determination of fy / fv- 



2 



1.2 How important are the couplings fv and /y? 

To exemplify the importance of having a good handle on fy/fv, we give the LCSR expres- 
sion for the form factor which dominates the B — > plv decay rate, namely [4] 



where M 2 is the so-called Borel parameter. The lower limit in the integral, uq, contains 
information about the energy of a threshold, So, above which quark-hadron duality is 
assumed. Using the input parameters of ref. [4], at q 2 = one has u$ ~ 0.65, whereas for 
q 2 = 10 GeV 2 it is u$ ~ 0.5. The objects under integrals are the light cone distribution 
amplitudes, which are the functions of u, the fraction of the longitudinal momentum of 
the meson carried by one of the meson's valence quarks. Indices T and || denote the 
polarization states of the p-meson. In the asymptotic limit, /i — > oo, these functions are 
^asymp.^ _ _ ^ [3J QCD corrections to (/>!^jj mp ' (w) are conveniently accounted for by 
expanding these functions in the basis of Gegenbauer polynomials and then by computing 
the first few moments. In order to make the above sum rule (1) as simple as possible we 
did not include radiative nor higher twist corrections. They can be found in ref. [5]. 

The first term in eq. (1), which is numerically dominant, is proportional to /J / f p . To 
show the importance of that ratio we first set it to one, as has been done in [3], and then 
to ' I f p ~ 0.75, as typically found in recent QCD sum rule analyzes [8, 9]. We plug those 
values in eq. (1), make the ratio of the two cases and plot it in fig. 1. From that plot and 
by varying the parameters as in ref. [4], we see that the form factor Ai(q 2 ) gets shifted by 
15-20% . Therefore a reliable estimate of fj/fp, as to reduce the overall uncertainty on 
Ai(q 2 ), is important. 

2 Definitions. Lattice details. Direct results 
2.1 Vector meson couplings 

For a generic vector meson V, with valence quark content (q'q), the couplings fv and fy(fJ>) 
are defined through the following matrix elements 



where p and are the momentum and the polarization vectors. These definitions refer 
to Minkowski space where we take = (i/2)[y tl , j v ]. Notice that since the anomalous 
dimension of the tensor current is different from zero, the coupling /y(//) depends on the 
scale at which the corresponding current is renormalized. 




(1) 



(0|g(0)7V(0)|V(p,A)) = /ymye£, 



(2) 
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2.2 Relevant correlation functions 



To estimate the values of the coupling fy(fj) for various light vector mesons V, we will 
compute its ratio with fy. This is beneficial because several systematic uncertainties, such 
as the dependence on the lattice spacing, are likely to cancel to a large extent. To compute 
fvilAI fv h is sufficient to consider the vector meson at rest (p = 0) and the following 
ratio of correlation functions in Euclidean space ( "£"' in the following) 

X 

where on the r.h.s. we chose the second operator to be the vector current. Any other 
local operator with the quantum numbers J PC = 1 would also be a good choice. The 
hat symbol indicates that the corresponding current has been improved and renormalized, 
which we discuss in the next subsection. By inserting a complete set of states and go- 
ing towards large time separations between the source operators, the lowest meson state 
dominates and we have 



(4) 



A 



A 



In eq. (4), we have Euclidean operators, whereas the couplings, as defined in eq. (2), refer 
to Minkowski space ( "M" ) . Using the fact that 

T f = (iq(x)a 0i q(x)) E - (q(x)a°*q(x)) M , Vf = {q{x) liq {x)) E - -i(q(x)Yq(x)) M , (5) 
it becomes clear that the plateau in the ratio (3), (i.e. eq. (4)) gives R(/i) = fvil 1 )/ fv- 



2.3 Improvement and renormalization 



Our lattice study is made by using non-perturbatively 0(a) improved Wilson fermions, 
where a stands for the lattice spacing. The improvement of the operators considered in 
this letter, T 0i (/j,) and V i: is made through [10] 



fi(n) = 4 0) M(1 + b T am q ) 



T£ (a) + c T (a)<9o^(a) 



Vf = Z ( y\a)(l + bv am q ) 



V»(a)-c v (a)d T E t (a) 



(6) 
(7) 



c r,v(a) ensure that bare currents computed on the lattice are free of 0(a) lattice artifacts; 
the constants Zy\a) and Z^\p,a) provide the matching to the continuum local opera- 
tors and also their renormalization in the chiral limit, whereas the quark mass dependent 
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artifacts of 0(am q ) are subtracted by setting the coefficients by (a) and bx{a) to their val- 
ues determined non-perturbatively and by using boosted perturbation theory, respectively. 
The numerical values of all of the above constants are listed in table 1, where we also give 
the main features of our lattices 1 . 



q _ a 1 „2 

P — o/fi'o 


o.u 


fi 9 


4 


CSW [1-Jj 


1 7fiQ 

i. ( oy 


1 /I 


1 .OZO 


/~ 3 x T 

±j A ± 


Ifi 3 v ^9 

-LU A 


94 3 y fi4 

A Ut: 


O^J A l U 


-44- cc\w T 


^oo 

O WW 


900 

Zj w w 


1 ^0 

-L tJW 


K -i 


n 1 33^ 


1 33Q 


1 347 


^2 


W . X ooo 


n 1 S44 


n i 34 q 


^3 


W . -LOtiW 


n 1 ^4Q 


W . -LOO -L 








1 






1 ^4 


1097 


ct [15] 


0.07 


0.06 


0.05 


cy [15] 


-0.11 


-0.09 


-0.08 


V [16] 


1.47 


1.41 


1.36 


b T [17] 


1.23 


1.22 


1.21 


4%) [12] 


0.766(2) 


0.775(2) 


0.795(3) 


4 0) (A*a = l) [12] 


0.833(2) 


0.847(3) 


0.867(6) 



Table 1: Summary of our lattice details. We also give the values of the improvement coefficients (with 
the corresponding references), and the values of the lattice spacing relative to the parameter r « 0.5 fm. 

2.4 Direct lattice results 

Information about our lattices is provided in table 1. Our data is obtained at three lattice 
spacings. The improvement ensures that the effects linear in lattice spacing are not present 

1 Please note that our lattice data refer to the same runs as those used in refs. [11, 12], where additional 
information can be found. 
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and therefore the extrapolation to the continuum limit (a — > 0) is expected to be smoother. 
In table 2 we collect results directly obtained for our lattices. 







Hi 


«2 


«3 


«4 


3 = 6.0 


rap/my 


0.743(6) 


0.709(8) 


0.682(10) 


0.650(12) 




fi/fv 


0.739(5) 


0.733(7) 


0.730(10) 


0.728(14) 




mpjmy 


0.807(4) 


0.751(7) 


0.661(11) 


0.575(16) 




fi/fv 


0.761(4) 


0.744(5) 


0.731(9) 


0.736(15) 


3 = 6.4 


mpjmy 


0.750(8) 


0.707(9) 


0.648(11) 


0.563(15) 




fi/fv 


0.747(4) 


0.737(5) 


0.728(7) 


0.723(15) 



Table 2: Ratio of the masses of the pseudoscalar and vector mesons consisting of two degenerate quarks 
of mass corresponding to m specified in table 1. We also give the values of the fy/fv ratio for each m, 
where the renormalization scale is set to /j = (l/a) p. 

In forming the ratio (3) we used time reversal to symmetrise the correlation functions 

as 

C T v(t) - \[C T v{t)-C TV {T-t)\ , 

Cyy(t) - l - [Cyy(t) + Cyy(T - t)} , (8) 

where t — 0, 1, . . . T/2. On the plateaus, the symmetrised correlator produces the usual 
tanh-factor, and thus the final fit form is 

f T 

R(t; l/a) = ^ tanh [m v (T/2 - *)] . (9) 
jv 

In improving the currents (see eq. (6)) we use the symmetric discrete derivative i.e. dog(t) = 
[g(t + 1) — g(t — l)]/2, where g(t) is the generic Green function computed on the lattice. 

Finally the fit is made using the same intervals with which the vector meson masses 
have been extracted (see ref. [11]), namely 

= 6.0 -> te [11,23] , 
P = 6.2 -> te [12,28] , 

= 6.4 -> te [17,28] . (10) 

The results are listed in table 2, and an illustration of the ratios R(t; l/a) is provided in 
fig. 2. 
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Figure 2: Ratio R'{t; 1/a) = fl(t; l/a)/tanh [my (T/2 - *)] (see eqs. (3,9)) which on the plateau leads 
to fy/fv- Fit to a constant is also shown. We show the signals for all three values of the lattice spacing, 
and we choose k = ki specified in table 1 for each (3. 
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3 Physical results 



3.1 Running 

We first need to run all the ratios from the scale fx = 1/a to the conventional fx = 2 GeV. 
In our computation we used the constant Zx(fx) determined non-perturbatively in the 
(Landau)RI/MOM renormalization scheme [12] which, up to NLO in perturbation theory, 
is the same scheme as MS(NDR). The scale dependence of the tensor bilinear, and thus of 
Z T (fx), is obtained by solving the renormalization group equation at NLO: 



-It Mi*)) , 



dlog/x 

where the anomalous dimension at NLO in perturbation theory reads [18] 



(11) 



It Ml*)) 



<o) <xM (1) 
4vr +7t 



An 



*°M + A (5 43 - 26n F ) 
3 4vr 27 V ; 



(12) 



At the same order, the running coupling is governed by 



da s (fx) 



-2f3 a s (fx) - ^-a 2 s (fx) , 
a log fx 2n 



(13) 



with the beta function coefficients being /3o — 11 — 2n F /3, and f3\ = 102 — 38np/3. Since 
fy(fx) inherits the scale dependence from the tensor operator, it is completely equivalent 
to solve eq. (11) in terms of fy(fx) instead of the matrix element (T^ (fx)). We finally have 



a»(/io)/ 



= KM 



X 



1 + Jl 



ajyfx) - a s (fi ) 

47T 



(14) 



where 



Jt = 



2 12411 - l2mn F + h2n 2 F 
9 (33 - 2n F ) 2 



Since we work in the quenched approximation, we set n F = 0, and we run to fx — 2 GeV 
by using A^l =0 = 0.250(25) GeV. The values of the lattice spacing are computed by using 
the conservative estimate a^ 6 = 2.0(1) GeV. The other a^ 1 are related to a^ 6 by using 
the ratios of a/r [14]. The net effect is that 



W[2 GeV, a' 1 ] = 1.000(2)(0) 



,1.010(2)(1) ,1.019(2)(1) 

,8=6.0 13=6.2 



13=6.4 



(15) 



where the first error comes from the 5% variation of ap= 60 , and the second one from 10% 
uncertainty on A^l -0 . As we can see those uncertainties are a few per mil and are thus 
negligible. 
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3.2 /J(2 GeV)//„ /£.(2 GeV)//** and fj(2 GeV)//, 

To get the physically relevant results (corresponding to V — p, if*, 0), we fit our results 
to the form 

#(2 GeV) _ , 2 , _ , , 4 



/ 



= a + ai(omp) +a 2 (o^p) , (16) 



where (amp) stands for the pseudoscalar meson mass that is computed by using the same 
quarks (non-degenerate in mass) as the ones used to compute the ratio on the l.h.s. (see 
table (2)). As our central values, we chose to quote the results obtained through the 
linear fit («2 = 0), whereas the quadratic fit is used to assess systematic uncertainties. 
The physical results, corresponding to V = p, K*, 0, are obtained by choosing (amp) 2 = 
(am^) 2 , (am K ) 2 , [2(am K ) 2 — (am T ) 2 ], respectively. The values of (am w ) 2 and (am K ) 2 are 
obtained from the fit of our data to 

(am v ) = 5 + 5i(am P ) 2 , (17) 

and the intersection of this form with the curve (amy) = r\J (amp) 2 , where r = 5.6, 1.8, 
i.e. (m p /m 7r ) phys ., and (m K * /m^phys., respectively. This extrapolation-interpolation to 
physical mesons is shown in fig. 3 and the results are presented in table 3. Our final results 
are obtained through a linear extrapolation to the continuum limit (a — > 0), i.e. 

tfecvn fm^n) (0)+B(a/ro? , (18) 



fv ! V / 

where on the l.h.s. we use our results as obtained at three values of the lattice spacing, the 
values of (a/r ) [14] on the r.h.s. are given in table 1, and B is a fit parameter. Notice that 
since our calculation is improved at 0(a), we omit the term linear in a when extrapolating 
to the continuum limit. The continuum extrapolation is illustrated in fig. 4 and our final 
results are 

fj(2 GeV)/f p = 0.720(24) (tj 6 ) , 
/£.(2GeV)//*. = 0.739(17) , (19) 
/J(2 GeV)//* = 0.759(9)(0). 



3.3 A brief discussion of the systematic uncertainties 

Many systematic uncertainties cancel in the ratio of the vacuum-to-meson couplings com- 
puted here. For example, the uncertainty due to the lattice spacing is only present in the 
evolution of the tensor current from the renormalization scale fj, — 1/a, down to /x = 2 GeV. 
As we showed, that uncertainty is at the level of a few per mil, thus negligible. 

Our central results at each value of /3, presented in table 3, are obtained by fitting 
our data to the form (16) linearly (a 2 = 0). The second errors reflect the difference 
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Figure 3: Ratios fy(2 GeV)/fy that are directly accessed in our lattice study (denoted by squares) 
are fit linearly (solid line) and quadratically (dotted line) according to eq. (16). To make the abscissa 
dimensionless we defined x = (amp) 2 / (am ss ) 2 , where (am ss ) 2 = 2(am^) 2 — (am^) 2 . Resulting ratios for 
the p, K* and <p mesons are marked by circles. The ratios obtained from the quadratic fit (gray circles) 
are slightly shifted to the left as to make them discernible from the ones obtained through the linear fit 
(red circles). From left to right they correspond to p, K* and <j> mesons. 
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Figure 4: Extrapolation to the continuum limit (a — > 0). 
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(a/r ) 2 


/J(2 GeV)// p 


/£.(2 GeV)//*. 


/J (2 GeV)/% 


0.035 


0.708(23) (±Sg) 


0.720(16) 


0.731(10)(0) 


0.018 


0.717(14) (t* 5 ) 


0.731(11) 


0.744(8)(0) 


0.011 


0.715(17) 


0.733(11) (lg) 


0.750(6)(0) 



Table 3: The values of the ratio of the vacuum-to-meson couplings mediated by the tensor vs. vector 
current, renormalized at /i = 2 GeV in the (Landau)RI/MOM scheme, which at NLO in perturbation 
theory is the same as the MS(NDR) renormalization scheme. 



between those and the values obtained through the quadratic fit (16). These errors are 
quite significant for the p-meson because that meson is reached through an extrapolation. 
On the other hand, they are small for the K*- and 0-mesons. 

The systematic errors induced by the presence of the chiral logarithms in extrapolating 
to the chiral limit are unlikely to affect our ratios because the same hadron state appears 
in both the numerator and the denominator. One may also speculate that for the same 
reason quenching effects should be small, but we prefer to consider our results as quenched 
and leave the unquenching issue to future lattice studies. 

In the continuum extrapolation we use a linear fit in the square of the lattice spacing 
since our results do not suffer from the 0(a) artefacts. In fact all the renormalisation 
and improvement constants are known non-perturbatively, except for 6 T which is obtained 
from 1-loop boosted perturbation theory. Since the residual 0(a)-effect is oc am q a 2 s and 
we work with the light quarks ((am,) < 0.07 [11]), this can be safely ignored. 

4 Tensor coupling and the comparison with other the- 
oretical predictions 

To get the absolute values of the tensor couplings we will now use our ratios from eq. (19), 
and multiply them by the physical values for f p ,K*,4> extracted from experiments, i.e. 

(rp . \ latt. 

! ' j "f ' a-,* (20) 
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4.1 Our results for ! pK *^ GeV) and comparison with other lat- 
tice calculations 

The experimental values for the charged vector meson couplings f p ±, f K *± can be extracted 
from the data for r-lepton decays [21]: 

BR(t -> p-v T ) = 25.0(3)% , BR(r -> K*~v T ) = 1.29(5)% . (21) 

By using the tree level expression 



BR(T _ y- Vr) _ SfcgfcE*** A *V l + ^V, (22) 



'V 



and = 0.9735, \V US \ = 0.220, r T = 290.6 ps, m T = 1.777 GeV, m K * = 0.892 GeV, and 
G F = 1.1664 • 10~ 5 GeV" 2 , all taken from ref. [21], we obtain 

f°T ~ 208 MeV , f^ ± -217 MeV , (23) 

where we display only the central values. Notice that the value of f p ± is consistent with 
f p o = 216(5) MeV, obtained from e + e~-annihilation. 

The constant can be obtained from the width T l po^ e + e - = 1.32(4) keV, which we 
combine with a" 1 = 137.036, = 1.019 GeV [21] and get 2 



/r=fr?^r,u e+c - ) .. 233 MeV. (21) 



1/2 

-I' „ ^ \ 

We then insert the above values for j e p x ^* $ m (20) and arrive at our estimate for the 
tensor couplings: 

/J±(2GeV) = 150(5) MeV = 152(7) MeV , 

/|,±(2GeV) = 160(4) (tj) MeV = 161(4) MeV , (25) 

/J(2 GeV) = 177(2) (0) MeV = 177(2) MeV . 

The result for fj (2 GeV) agrees also very well with the estimate of ref. [7] in which 
this constant has been computed on the lattice, at three values of the lattice spacing, in 
the quenched approximation. In that study the perturbatively evaluated renormalisation 
constant has been used. Their new calculation [19] includes also the non-perturbative 
determination of Zt((j), and after extrapolating to the continuum limit they quote 

fj (2 GeV) = 150(4) MeV . (26) 

Although the systematic uncertainty of that result is still preliminary, we note a very 
pleasant agreement with our value given in eq. (25) 3 . 

2 In the previous version of this paper, eqs. (22) and (24) contained typos. Numerical results however 
were correct. We thank Patricia Ball for drawing our attention to this. 

3 We are indebted to Gerrit Schierholz for pointing out the ref. [7] to us and for communicating their 
new results. 
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4.2 Comparison with the QCD sum rule results 



The QCD sum rule (QSR) estimates [5, 8] are made at the renormalization scale \x = 1 GeV. 
We run those values to fi = 2 GeV, by using eq. (14) and A^ =3 = 338(40) MeV [20], i.e. 
ft (2 GeV) = 0.9667(6) x ft(l GeV). Thus the QSR estimates are 

ft {2 GeV) = 155 ± 10 MeV , 

ft, (2 GeV) = 179 ± 10 MeV , (27) 

ft (2 GeV) = 208 ± 15 MeV. 

Note also that the coupling ft(fx) has been re-examined in ref. [9], essentially confirming 
the results of ref. [8]. By comparing the lattice results (25) with the QSR values (27), we 
see that the agreement of the two methods for ft is very good, quite good for ft* and less 
good for ft. Investigation of the source for that discrepancy is beyond the scope of the 
present paper. It should, however, be stressed that the LCSR yield results for full QCD 
while our results are quenched and therefore a perfect match is anyway not expected. 

5 Summary 

In this letter we presented quenched lattice results for the ratio of the couplings of the 
light vector mesons to the tensor and to the vector currents. Those ratios enter the LCSR 
analyses of the phenomenologically important heavy — > light meson semileptonic decay 
form factors. From the results obtained with 0(a) improved Wilson quarks, and with high 
statistics data, we were able to extrapolate to the continuum limit. The resulting ratios 
are then multiplied by the experimentally measured vector meson couplings to arrive at 
the tensor couplings. Comparison with the prediction obtained from the QCD sum rule 
analysis indicates good agreement for ft K *{2 GeV), whereas the agreement is not good in 
the case of ft (2 GeV). Our result for ft K *{2 GeV) also agrees very well with the lattice 
estimate of refs. [7, 19]. 
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